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In a previous study {Comp. Phys. Com. (2005) 169, 139-143), we clarified the dependence 
of the phase structure on the hydrophilicity of an amphiphilic molecule by varying the interaction 
potential between the hydrophilic molecule and water (haw) in a dissipative particle dynamics 
(DPD) simulation using the Jury model. In the present paper, we perform another DPD simulation 
using the previous model to investigate the dependence of the interaction potential between adjacent 
hydrophilic groups on the phase structure. By varying the coefficient of the interaction potential 
^ between adjacent hydrophilic groups qaa (oaa ~ 15, 25, 40, and 250) at a dimensionless temperature 

'^'h I of T = 0.5 and a concentration of amphiphilic molecules in water oi (f> — 50%, hexagonal (aaa = 

fvj ■ 14, 25, 40) and micellar (aaa = 250) phases were observed. In comparison with the previous results, 

the dependence of the A-B dimer's shape on haa was determined to be weaker than that on aAw- 
Therefore, it is concluded that the solvent water W plays an important role in aggregation of the 
A-B dimers. 
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Amphiphilic molecules have many degrees of freedom in the structures they adopt. However, when the temperature 
and concentration of an aqueous solution of the amphiphilic molecules are fixed, the molecules are restricted to a 
Q I certain molecular shape and aggregate to form a variety of mesoscopic structures, for example, micellar, lamellar and 
Q ■ hexagonal phases. 

In previous papersP,l3lj we have investigated the amphiphilic molecule hexaethylene glycol dodecyl ether (C12E6), 

. ■ a popular surfactant in water that forms a variety of self-assembled structures. The phase structure of C12E6 was 

^ ' investigated by Mitchell^] in 1983. In our work, we clarified;^] the dependence of the phase structure on hydrophilicity 

\j^ by varying the interaction potential {a aw) between the hydrophilic molecules and water in a dissipative particle 

^\ [ dynamics (DPD) simulation using the Jury modelH|. 

04 I In the present paper, to determine the relative dominance of different microscopic interaction potentials in de- 

J^ ' termining the mesoscopic phase structure, we investigated the dependence of the phase structure on the interaction 
potential between adjacent hydrophilic groups caa- The details of the DPD simulation algorithm and models can 
not be included due to lack of space but are the same as reported previously^ with the exception of the interaction 
potentials oaw and caa- The differences with the previous simulation are described in Section 2. 
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II. SIMULATION METHOD 

5 ■ We used the same DPD model and algorithm jj,|3 as described previously[ll|2|, with a modification of the conser- 

fj vative force F^j between particles i and j, given in the present paper by 

^c ^ f a»j(l - nj)n,j if ry < 1, ,^. 

'■J 1 a rij > 1, ' 

where n^ = (r^ — rj)/\ri — rj\, and r^ is a position vector for particle i. Coefficients Uij in Eq. ^ denote the 
coupling constants between particles i and j. The numerical values of a^ are given by Tabled 

III. SIMULATION RESULTS AND DISCUSSION 

To demonstrate the dependence of the mesoscopic phase structure on the interaction potential between adjacent 
hydrophilic groups, we varied the A-A interaction potential coefficient such that oaa — 15,25,40, and 250. 
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FIG. 1: Structures formed for each potential coefficient, baa ~ 15(a), 25(6), 40(c), and 250(d). We set T — 0.5 and (f> = 50% 
during simulation. 



The aggregation structure for each value of aaa is shown in Fig. ^ To quantitatively classify the phase structure, 
we also plot the radial distribution function of the solute particles for each oaa in Fig. [3 The figures show that 
hexagonal phases (oaa = 15,25,40) and micellar (aaa = 250) phases are formed depending on the coefficient aaa of 
the A-A interaction potential. 

The distance R between a hydrophilic group in an A-B dimer and another group in the nearest A-B dimer (Fig. ^ 
is shown as the first peak R{aAA) of gAA{r) in Fig. |2Ia). On the other hand, the intramolecular distance between A 
and B, that is I in Fig. O is shown as the first peak in Fig. |2Ib). Wc plotted the dependence of the distance R and 
the length / on the interaction potential oaa in Fig. ^a). Figure QJa) indicates that the dependence of the length 
Z(aAA) on AAA is weaker than that for R. We also plotted the dependence of 1/{R^I) on the interaction potential caa 
inFig.Hb). 

Here we consider the packing parameter introduced by Israelachvilij^ |3, l^l- According to this description, the 
hexagonal and micellar phases correspond to 1/2 < p < 1/3 and p < 1/3, respectively. The packing parameter p is 
considered to be proportional to 1/{R^I). From Fig. Efb), p oc 1/{R?1) decreases as oaa increases. Therefore, it is 
found that A-B dimers form cylinders for small oaa values, and that AB dimers modify their shape from cylinders 
to cones by increasing the distance between the neighboring head groups as oaa increases. 

Lastly, we compare the present results with our previous workj^l, which investigated the dependence of the phase 
structure on oaw- In that previous work, we showed that the A-B dimer modifies its shape with changing hydrophilicity 
oaw ■ It is intuitively expected that the shape of the A-B dimer would depend more strongly on the parameter caa 
than the parameter cawj because oaa affects the distance between the hydrophilic groups more directly than oaw- 
However, contrary to our expectation, the present simulation results show that the dependence of the A-B dimer's 
shape on caa is weaker than that on caw- Based on this result, we conclude that the solvent water W plays an 
important role in determining the phase structure of the A-B dimers in aggregation. 
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FIG. 2: Radial distribution functions for the solute particles vs. distance between the two particles r for baa = 15, 25, 40 
and 250. (a) The A-A radial distribution function QAAir). The first peak R{aAA) of each curve corresponds to the A- A 
distance between two adjacent AB dimers. i?(15) = 0.889, 7?(25) = 0.889, i?(40) = 0.904, i?(250) = 0.963. (b) The A-B radial 
distribution function gABir). The first peaks correspond to the length of the A-B dimer I in Fig|3| Z(15) = 0.326, Z(25) = 
0.326, Z(40) = 0.326, Z(250) = 0.296. 
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TABLE I: Table of coefficients a,;j showing 
dependence on particle type for particles 
i and j, where W is a water particle, A 
is a hydrophilic particle, and B is a hy- 
drophobic particle. By varying the coeffi- 
cient AAA between adjacent particles A and 
A, aAA = 15,25,40 and 250, the depen- 
dence of the phase structure on the A-A 
interaction potential can be clarified. 




FIG. 3: Schematic diagram of two adjacent 
molecules. A gray ball and a twisting black 
line are used to denote the hydrophilic and 
hydrophobic parts, respectively, of an am- 
phiphilic molecule. The parameter R is the 
distance between the hydrophilic groups, 
and I is the "maximum effective length" of 
the hydrophobic tail. 
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FIG. 4: Shape parameters of two adjacent molecules vs. the A-A interaction potential coefficient oaa. (a) The A-A distance 
between two adjacent AB dimers R and the length of the hydrophobic tail I are shown, (b) The inverse of 7?^i, which is 
proportional to the packing parameter p, is plotted. 



IV. CONCLUSION 

We demonstrated the dependence of the A-B dinicr's aggregated structure in solvent water W by varying the A-A 
interaction potential caa in a dissipative particle dynamics simulation. The present simulation and the previous 
resultsj2| show that the OAA-dependence of the A-B dimer's shape is weaker than the aAw-dependence. Therefore, it 
is concluded that the solvent water W plays an important role in forming the phase structure of the A-B dimers in 
aggregation. 
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